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Two-loop corrections for the standard Abelian Nambu- Jona-Lasinio model are obtained with the 
Optimized Perturbation Theory (OPT) method. These contributions improve the usual mean-field 
and Hartree-Fock results by generating a 1/N C suppressed term, which only contributes at finite 
chemical potential. We take the zero temperature limit observing that, within the OPT, chiral 
symmetry is restored at a higher chemical potential fi, while the resulting equation of state is stiffer 
than the one obtained when mean-field is applied to the standard version of the model. In order 
to understand the physical nature of these finite N c contributions, we perform a numerical analy- 
sis to show that the OPT quantum corrections mimic effective repulsive vector-vector interaction 
contributions. We also derive a simple analytical approximation for the mass gap, accurate at the 
percent level, matching the mean-field approximation extended by an extra vector channel to OPT. 
For fi > /i c the effective vector coupling matching OPT is numerically close (for the Abelian model) 
to the Fierz-induced Hartree-Fock value G/ (2N C ), where G is the scalar coupling, and then increases 
with /i in a well-determined manner. 
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I. INTRODUCTION 



Effective models are extensively used to understand the physics of strong interactions. This is particularly true in 
the study of the structure of the phase diagram of Quantum Chromodynamics (QCD) , specially in the low temperature 
(energy) and high density region, which is supposedly almost inaccessible through the present day lattice techniques 
that make direct use of QCD. In this respect, the use of effective models for quark interactions, like for example the 
Nambu- Jona-Lasinio (NJL) type of models has proven to be extremely helpful to improve the understanding of 
the phase structure for strongly interacting matter. 

With NJL type of models the phase structure of QCD has been explored mostly in terms of the well known 1/N C 
expansion, whose leading contribution represents the large- N c (LN) approximation and also corresponds to the mean- 
field-approximation (MFA) [J. Recently, the alternative Optimized Perturbation Theory (OPT) method has been 
applied to the NJL model Q and results beyond the MFA have been explicitly obtained for quantities related to the 
QCD phase structure. This application has been extended to the strong coupling and small current mass regime in 
order to investigate the critical line on the chemical potential-current mass plane [4|. 

For example, it has been observed that, at T = 0, the OPT predicts that the first-order phase transition takes place 
at a chemical potential value that is higher than the one predicted by the MFA, thus, chiral symmetry breaking is 
enhanced. At this point it is interesting to remark that this trend is also observed when the MFA is applied to the 
NJL augmented by a repulsive vector- vector interaction term, — Gy {^Pj^ip) 2 - It is then plausible to imagine that the 
OPT 1/N C correction captures some of the physics associated with this kind of interaction, which plays a major role 
at finite densities. The aim of the present work is to address this possibility by performing a numerical comparison 
between the OPT results, obtained from the standard NJL model, and the MFA results, obtained with the extended 
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version of model, such that the role of the quantum corrections generated by the former approximation can be better 
understood. 

Although of academic interest only, the Abelian version of the NJL is well suited for such a task, since in this 
case, due to cancellation of terms between scalar and pseudo-scalar contributions at the exchange level [2|, only a /i 
dependent contribution survives, allowing for a clearer comparison. Therefore, the temperature does not play a very 
important role in the present study and we, consequently, can restrict the analysis to the T = case. By adjusting 
the value of the vector-vector coupling in the MFA we will perform a numerical analysis exhibiting that the OPT 
results can be approximately well reproduced, near criticality and quite beyond, at the value Gy ~ G/(2N C ), where 
G is the scalar coupling. As we shall discuss, this could be anticipated by comparing the OPT and MFA relations for 
the free energy density. We then argue that one of the OPT main effects is to induce a new type of physics within the 
standard version of the NJL, by radiatively generating a vector- vector type of 1/N C correction that mainly affects the 
finite fj, results. The possible generalization of these findings and their consequences for general four-fermion effective 
models, such as the non- Abelian NJL model as well as the Gross-Neveu model, and perhaps even for QCD, is currently 
being investigated. In this respect, the present investigation, that is carried out in this simple case, is important in 
its own right, since it offers the necessary framework to pursue such generalizations. We should also notice that the 
particular value Gy — G/(2N C ) at the level of our approximations just corresponds to the Fierz-induced Hartree-Fock 
(HF) approximation, which is not surprising since the OPT at first order involves corrections topologically similar to 
Fierz exchange terms. However, when considering the /i-dependence, as we shall illustrate the OPT incorporates more 
corrections than simply Fierz-exchange terms, and the Gy value best fitting near criticality is merely a numerical 
accident, only valid for the Abelian model. An additional bonus of our study is thus to produce a simple general 
analytical comparison of the MFA with vector-vector interaction approximation with the OPT one, from which we 
will be able to show the differences and advantages of the OPT at first-order over the HF approximation. 

The remaining of this work is organized as follows. In Sec. |TT] we derive the free energy density using the OPT 
formalism. In Sec. IIIII we perform a numerical comparison between the OPT and the MFA results for the thermo- 
dynamical quantities at T = 0. The analytical comparison of the OPT results with those coming from MFA with 
Gy ^ is studied in Sec. IIV1 Our conclusions are presented in Sec. |Vj 



II. THE OPT FREE ENERGY FOR THE ABELIAN NJL MODEL 



The simplest version of the Abelian NJL model is described by a Lagrangian density for fermionic fields given by 1 



C = i> (t? - m c ) ip + G [{^f + (^ 7 5^) 2 ] , (2.1) 

where ip in general represents a iV c -plet quark of just one flavor with current mass m c . Therefore, here, we consider 
the one-flavor model with (global) U(l) symmetry. In the OPT interpolation prescription Q applied in the case of 
four-fermion theories [3, HI, El, one starts by deforming the original theory with the replacements m c — > m c + (1 — S)i] 
and G — > SG in Eq. (|2.ip where S is the alternative perturbative expansion parameter that determines the order at 
which the OPT expansion is carried out and 77 is an arbitrary mass parameter to be set in a variational way through 
an appropriate optimization procedure at a given order in the OPT. The typical optimization procedure used in the 
many previous applications and that we also adopt here, is the Principle of Minimal Sensitivity (PMS), defined by Q 



dr] 



rj, 8=1 



(2.2) 



where V^ k ' is some physical quantity calculated up to order k in the OPT. From the above prescriptions to build 
the interpolated model and rewriting the quartic interaction in (|2. L|) through a Hubbard-Stratonovich transformation 
introducing auxiliary fields, a and 7r, we obtain the interpolated Lagrangian density in the OPT formalism, 



C = 4> {1$ - (m c + 77) + 5 {{r, -a)- tj 5 ir]}^ - 6^ (a 2 + tt 2 ) . (2.3) 

Note from the interpolated theory Eq. (|2.3I) that the original Lagrangian is recovered for 8 = 1, but at any finite 
order k in the OPT, any perturbative result evaluated with Eq. (j2.3p becomes dependent on 77. Fixing it through 
the optimization procedure (|2.2p . new kinds of contributions beyond MFA are embodied in the non-trivial optimal 
?7-dependence. 



3 



For the studies involving the thermodynamics and phase structure of a given field theory model, the most appropriate 
quantity to be optimized has been shown 5, 6] to be the effective potential (or free energy density), J 7 . Applying 
the optimization procedure (|2.2|) to T evaluated to some order k in the OPT, we can then access the system's 
phase structure in a nonperturbative way beyond MFA. From the interpolated Lagrangian density, Eq. (|2.3I) , we can 
immediately read the corresponding Fcynman rules and then perform a perturbative expansion of the free energy 
density in powers of 5. Then, the free energy density in the <j c direction (using that (a) = a c and (tt) = 0) reads [3j 
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All the integrals in Eq. (|2.4I) are to be interpreted in the Matsubara's finite temperature formalism with 



(2.4) 
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(2.5) 



and quadri-momenta given as p = (iu) n + /_t, p), where w„ = (2n + 1)ttT, n = 0, ±1,±2, are the Matsubara 
frequencies for fermions. By performing the sum over these frequencies (see e.g. Ref. [3j), one obtains 
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where we have defined, for convenience, the following basic integrals: 
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(2.9) 



where E 2 = p 2 + (77 + m c ) 2 . Here, as in Ref. only the divergent integrals occurring in Eqs. (|2.7|) and (|2.8j) are 
regularized by a sharp non-covariant three-dimensional momentum cutoff, A, so that the Stefan-Boltzmann limit is 
observed. Note that the second term in Eq. (|2.4I) corresponds to a gas of free fermions, whose mass has been dressed 
by rj, while the third term represents tadpole type of contributions, proportional to the quark condensate, The 
last term in Eq. (|2.4[) . which is proportional to the square of the quark number density, (■0 + '0) 2 , comes from the 
two- loop terms in a next-to- leading order correction in the 1/N C expansion, which can be easily seen by redefining 
G — > G/N c . This term is responsible for generating the vector type of corrections. 

In order to perform our evaluations, we need to consider the general PMS equation (|2.2p . which can be conveniently 
expressed in the form Q 



1 + (77 + 77l c ) 



drj 



2Gh—I 3 
dr\ 



0. 



(2.10) 



One observes from Eqs. (|2.6p and (|2.9[) that the first 1/N C correction 1 (represented by the last term) is only relevant 
when /i 7^ 0, as it vanishes for (i = 0. In the latter case, the PMS condition gives fj — a Cl merely reproducing the well 
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known MFA results. Thus, the relevant contribution arises when p ^ 0. Since we are mainly interested in the phase 
structure of the model and on the effects of the OPT on it, one basic quantity of interest is thus the thermodynamical 
potential, f2, whose relation to the free energy is given by = J-(cr c ). The optimized order parameter, a Cl is 
determined from the gap equation generated by minimizing JF with respect to the classical field, o~ c . From Eq. (|2.6[) 
we obtain that 

a c = AGN c (fj + m c )I 2 (p,T,fj) , (2.11) 
which can be inserted in the PMS equation to determine the optimized pressure. 

III. MATCHING THE OPT CORRECTIONS IN TERMS OF A VECTOR INTERACTION 

The effects due to the OPT first-order corrections for cold and dense quark matter can now be studied for the case 
of finite density effects (p ^ 0) at zero temperature which is easily achieved by taking the limit T — > in the thermal 
integrals [1]. Physically, this situation is relevant in studies related to neutron stars for example. This is also the 
regime where lattice techniques face more problems with the sign problem. Having obtained the optimized pressure 
one can easily obtain the quark density number, p = dP/ dp, and by taking into account both the gap and PMS 
equations. The energy density is given as usual by e = —P + pp and allows us to obtain the equation of state (EoS). 

Now, in order to carry out comparisons with the MFA let us recall that, as emphasized in Refs. ja-UjJ, the 
introduction of a repulsive vector- vector interaction term of the form -Gv^"^) 2 in Eq. (|2.1[) is also allowed by 
chiral symmetry. Such a term can become important at finite densities, generating a saturation mechanism, depending 
on the vector coupling strength, that provides better matter stability. Also, as far as phase transitions are concerned, 
it can influence the size of the first-order transition region and, hence, the location of the (tri)critical point. As already 
mentioned, the main goal of the present work is to compare the OPT results obtained at Gv = 0, with the MFA 
results for Gv i= and whose free energy density is given by |j[ 



where Mmfa is the mass gap in MFA. In this case, one has to simultaneously solve the equations 

P = p- 4G v N c I 3 (ii, T, Mmfa) , (3.2) 

and 

Mmfa = m c + 4GN c M M FAh (A, T, M MFA ) . (3.3) 

For our purpose, it is crucial to notice that the combined effects of p, defined by Eq. Q3.2p . and of the — 4Gy_/V 2 7| 
term in Eq. fl3.1|) . is to produce a net positive contribution approximately ~ 4GvN 2 i| (for p ^> M) to the the free 
energy [jjl, Il3|. This last term turns out to be of the same form as the OPT finite N c correction term, 2GN C 2|, so 
that one may guess that results obtained with this approximation should be close to those obtained with the MFA at 
Gy values of the order G/(2N C ), but not the same (see detailed discussion in Sec. II VI) . This is precisely what we will 
check now in a numerical fashion by considering the EoS. 

Following the standard procedure for the NJL model, we set the OPT parameter values to A = 620 MeV, GA 2 — 4.44 
with m c = 5.0 MeV, so that at T = one matches by convention th e "p hysical" values of f n — 92.4 MeV, = 
135 MeV, while predicting M = 341 MeV, and ~{4>i>) 1/3 = 191 MeV 0- One subtlety is that the m c value thus 
determined for OPT is strictly speaking slightly different from the corresponding MFA m c , when obtained consistently 
from fitting the same pion mass value: this is due to the fact that the OPT gives G 2 correction to the relation between 
m c and Q. More precisely, one finds ■m\ iFA ~ 5.36 MeV instead of m^ PT ~ 5.0 MeV for the same A = 620 MeV 
input. (In contrast the G value being obtained from f n is not changed from MFA to OPT, since for the U(l) model, 
fn receives no OPT corrections at G 2 order). Within the MFA, we will also use the following representative values 
for the vector coupling, G v = 0, 3G/(87V C ), G/{2N C ) and 3G/(4A^ C ). 

Figure [T] shows the quark effective mass as function of the chemical potential /i, obtained with the OPT and the 
MFA for the four different Gv values. First, one notes that the p c value at which chiral symmetry is restored, via 
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FIG. 1: Effective quark mass (at T — 0) as a function of the chemical potential. 



TABLE I: The critical chemical potential fi c (T = 0) in the OPT and in the MFA (for the different values of G\ 



Gv 





3G/(8iV c ) 


G/(2N C ) 


3G/(4N C ) 


OPT (G v = 0) 


(MeV) 


350.9 


356.3 


358.3 


361.9 


358.5 



a first order transition, increases with Gy in the MFA. The critical chemical potential value predicted by the MFA 
with Gy = G/(2N C ) agrees the better with the OPT prediction, as we can see also from the results shown in Tab. fl] 
The pressure subtracted by the ii — value, P(/i) — P(0), as a function of the chemical potential is shown in Fig. [5] 
We also see from the results that the OPT prediction is very close to the MFA with Gy = Gj (2N C ) near criticality, but 
the OPT curve is less steep than the MFA one and tends to deviate towards higher values of Gv as fi increases. Thus, 
one may expect that the baryon density predicted by the OPT to be smaller at high /i. This expectation is confirmed 
by the results shown in Fig. [3l where we show the baryon density as a function of the chemical potential. The same 
behavior is also reflected in Fig. HJ where is shown the energy density as a function of fi. Again, one sees that the 
OPT and the MFA with Gv = G/(2N C ) are in best agreement around criticality, as shown by the results in Table [H 
while at higher values of // the OPT result interpolates between this and higher effective Gv values. This behavior of 
the OPT results at larger chemical potential can be predicted if we compare the OPT expression for the free energy 
density, Eq. (|2.6I) . with the one in the MFA approximation, Eq. (I3.1[) . Besides of the nontrivial dependence of the 
expression, in the OPT case, with the variational parameter 77, there is also the additional contribution with the I2 
term in Eq. (|2.6p . All these contributions from the OPT tend to favor a comparatively larger value of Gv in the 
MFA case, for higher values of the chemical potential, when trying to match the OPT results with the ones within the 
MFA with Gv i= 0. In the next section we will trace this behaviour by deriving a more analytical connection between 
the OPT and such an effective mean field vector-vector coupling. Finally, this behavior of the OPT compared with 
the MFA case with a non vanishing value of Gy, also reflects on the EoS. The EoS results for the different cases are 
shown in Fig. [5] and serve to illustrate how the OPT corrections produces a stiffer equation of state for this version 
of the model 3. 



IV. A SIMPLE ANALYTICAL APPROXIMATION 



We shall try to get a deeper insight about our previous numerical results by examining a simple well-motivated 
approximation. As we will see, it will essentially explain the bulk of the OPT results versus the MFA ones with extra 
vector coupling, explaining in particular the rather accidental best OPT matching by the HF Gy — G/(2N C ) value near 
criticality. First, as it is well-known, the large- N c (or MFA) is also equivalent to the traditional Hartree approximation, 
which only incorporates effects from direct (tadpole) contributions, while exchange terms are also considered within 
HF. Indeed, the value Gy = G/(2N C ) is just the vector coupling one obtains by a Fierz transformation of the original 
interaction, Eq. (|2.1[) . Therefore, Eq. (|3.1[) with this value of Gy just corresponds to the HF approximation. Similarly, 
the two-loop induced first order OPT corrections to the free energy incorporate the structure given by 1 /N c exchange- 
graph of the Fierz type. Indeed for fj, — (and T = 0) the OPT results are generically consistent @ with HF ones 
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FIG. 2: The pressure (subtracted by the value at ii — 0) as a function of the chemical potential. 
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FIG. 3: Baryon number density, ps = p/3, in units of nuclear matter density (po = 0.17 fm 3 ), as a function of the chemical 
potential. 



(up to higher 1/N% corrections originating from the optimization of the mass rf). Note however that in the Abelian 
case, Mhf is identical to Mmfa, due to the exact cancellation of the 1/N C scalar and pseudo-scalar contributions. 
However, as far as the ^-dependence is concerned, OPT corrections also induce extra terms, as is already clear from 
Eq. (|2.10j) . Comparing the set of equations above one sees that nonperturbative information concerning a vector type 
of interaction is being taken into account by both approximations in a different way. This becomes clear if one recalls 
that the I3 integral roughly represents the quark number density, while the integral I2 roughly represents the 

scalar condensate, (ipip) ■ The last term of Eq. (|2 . 10[) explicitly displays that J3 is also considered in a nonperturbative 
fashion. 

In order to understand further beyond purely numerical results the origin of these differences, let us have a deeper 
look at the expressions for the MFA and OPT mass gap expressions, respectively given by Eq. (|3.3[) with (I3.2[) for 
the former (where I2 and -Z3 are to be evaluated at p), and Eqs. (|2 . 1 1 1) and p.lfjp for the latter. For fi > /i c a simple 
but crucial observation is that the behavior of both expressions is essentially driven by the explicit chiral symmetry 
breaking mass m c . Indeed, for m c = 0, one would have Mmfa{Gv)(^ > Ac) = and Mopt{^ > Mc) = consistently. 
Now, since m c is small as compared to all other scales, we may expect that to a good approximation the mass gap 
could be well approximated by a first-order expansion in M, next to be solved by a self-consistent equation (becoming 
in that case a trivial first-order equation): 



M ~ m c + 0(n, A, G, • • • )M + 0(M 3 ) , (4.1) 
where /?(/i, A, G, ■ ■ ■ ) is the coefficient obtained from this first-order expansion. We also note that the next-order term 
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FIG. 4: Energy density as a function of the chemical potential. 
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FIG. 5: Equation of state for quark matter (pressure as a function of energy density). 



in the expansion in Eq. (|4.1j) is of order M 3 . 
Explicitly, we obtain after some simple algebra: 



where from Eq. 



and 



N C G 



M M fa(Gv) ^m c + (A 2 - /2 2 ) M MFA + 0(M 3 ) , 



a = „ - 2 ^A 3 + o(am 2 ) 



(A 2 -/i 2 ) + 



3iV c (A 2 - M 2 ) 

respectively. The gap equation solutions are thus trivially given as 



M, 



OPT 



0(M 3 
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and similarly for Mmfa(Gv). Note that those approximations both turn out to be accurate at the percent level, as 
compared to the exact numerical solutions, for a large range of \i c < fi <C A values. This is largely due to the very 
suppression by higher 1/(A 2 — /z 2 ) powers of higher order coefficients of these expansions in powers of M . 

Now defining generically Gy = a G/N Cl a will depend essentially on fj, and this is what we would like to determine 
next. Before that, we observe that, at least for rather moderate values of Gy ~ 0(G), Eq. (|3.2p when solved exactly 
actually gives a moderate shift of p, with respect to [i, giving // slightly smaller than [i by only a few percent. So 
fi c is well approximated by restricting Eq. (|4.3[) at first iteration, i.e., first Gy order. This allows to derive a simple 
relation to determine the relevant coefficient a above for the Gy /G ratio, as a function of y, and the other parameters. 
From the matching Mm fa = Mopt, then from Eqs. (|4.2|) and (|4.4[) and using Eq. (I4.3[) at first iteration, but solving 
exactly the (quadratic) equation for a, after a simple algebra one finds 



which gives an explicit relation for the appropriate Gy values needed to match the OPT results for arbitrary y (only 
valid of course for /i > /i c , and also for /j < A in principle, since reaching too close to the natural cutoff does not 
make much sense physically within the NJL effective model). 

In deriving Eq. (|4.6[) it is assumed that both the OPT and MFA masses involve the very same current mass m c , 
which is thus valid if m c is a common input to both approximations. However, as discussed above in Section IlIIl 
the OPT m c value as determined consistently from fitting is slightly smaller than the corresponding MFA value. 
Although a rather modest change in m c it affects the results, due to the high sensitivity to m c for // > // c , making 
the OPT mass gap lower by the same amount according to Eq. (|4. 5|) above, as compared with the MFA mass gap for 
a given Gy value. This partly delays the increase in the effective Gy for increasing [i from Eq. (|4.6p . which can be 
easily taken into account and combined with the latter to extract a precise relation for Gy /G appropriately modifying 
Eq. (|4.6[) . Since m^ PT /m^ IFA ~ 0.93, this gives a moderate modification to the relation (j4.6[) in practice, making 
the corresponding effective Gy slightly smaller by about 10% for a given /i. This fits remarkably well our numerical 
results for the mass gap obtained with the exact expressions, with only a few percent errors. For /x ~ fj, c ~ 360 MeV, 
i.e. very near criticality, one finds GyN c /G ~ 0.53, which explains a posteriori why the results remain there very close 
to the HF value GyN c /G = 1/2. However, it is also clear from the previous derivation that this is to some extent 
a numerical accident, resulting from a rather fortuitous combination of the Abelian model factors, GA 2 ~ 4.44 and 
/i c /A values such that the term in parenthesis in Eq. (|4.6p is ~ 0.78, plus the small correction from mJ? PT /m^' IFA as 
above explained. Indeed for higher fi Gy jG grows first moderately and then more rapidly when [i increases towards A. 
For example, we obtain Gy(/i ~ 400 MeV) - 0.6G/N C and Gy{fi ~ 500 MeV) - 1.13G/iV c . This also illustrates the 
differences between OPT and HF or any other fixed Gy value, since OPT mimics an effective Gy(/i) value. Finally, 
the behaviour of the other physical quantities in Figs. 2-5 may be understood semi-analytically along the same line 
of reasoning, though we skip the details here. The pressure in Fig. 2 follows roughly the same trend as the mass 
gap. As for the other quantities in Figs. 3-5, the more pronouced increase of the effective OPT-matching Gy value 
at increasing /i values can be understood essentially by the enhanced /^-dependence from the definition of the density. 



In order to investigate the physical origin of the main effects produced by the OPT nonperturbative approximation, 
we have evaluated the free energy density for the Abelian NJL model. This version of the model is certainly less 
realistic than the non- Abelian ones, but being much simpler, it suits our purpose by providing an easier to analyze free 
energy. As usual, in the large- N c limit, the OPT reproduces the MFA result exactly, while at the next-to-leading 1/N C 
order as induced by OPT, the contributions from the scalar and pseudo-scalar channels cancel each other whenever 
H = (for this version of the model 0]). Nevertheless, the finite N c contributions may still be important for situations 
such as T = 0, fi ^ (dense cold quark matter) and T ^ 0, [i ^ (hot and dense quark matter). 

Considering the MFA predictions for the same version of the model one observes that the OPT prediction for the 
critical chemical potential value (/i c ), at T = 0, is higher. However, as it is well known, the MFA applied to the 
NJL in the presence of an extra repulsive vector- vector interaction (proportional to the coupling Gy) also predicts 
higher /i c values, as compared to the Gy = case. Then, comparing the mathematical structure of the free energy 
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density provided by the OPT, with Gy = 0, and the MFA, with Gy ^ 0, we show that the former approach radiatively 
generates an effective vector-like type of contribution which is 1/N C suppressed and, hence, does not appear in a large- 
N c type of calculation. This comparison hinted to the value Gy = G/(2N C ) as the one for which both approximations 
could furnish similar results in a reasonable range of fi values near criticality. This has been successfully verified in 
both numerical and analytical manner, showing that the OPT predictions, obtained with one less parameter (Gy), 
point out in the same direction as the ones produced by the MFA applied to the NJL with Gv ^ and, thus, with 
a larger parameter space. In particular, the agreement for the value of Gy — G/(2N C ) is best near criticality, but we 
have shown by a simple analytical approximation how OPT incorporates corrections for higher fi matching higher Gy 
values. Indeed it should be remarked that the very proximity, near criticality, of the OPT results with the simpler HF 
approximation, is not much a generic feature but a rather fortuitous accident of the very simple Abelian model here 
considered. In the SU (2) version of the NJL model, the OPT exhibits Q more departure from the HF approximation, 
and in other models also it generically captures corrections beyond HF (for instance in Ref. [l~5j results obtained from 
OPT at first order for the Gross- Neveu mass gap are very close to the full next-to- leading 1/N corrections). Indeed an 
interesting follow up beyond the scope of the present work could be to examine similarly the correspondence between 
a generalized OPT-NJL SU(N) and extra four-fermion interactions in MFA. 

The present application suggests that the OPT is able to bring new physics features, without introducing new 
parameters. This conclusion certainly deserves further investigation and we are carrying out a similar analysis for the 
case of general four-fermion theories in an attempt to estimate the consequences beyond the MFA for the QCD phase 
diagram with a special concern about the location of the critical point. These results will be presented elsewhere. 
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